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Abstract: Dirac fermions in the fundamental representation of SU(N) live on the surface 
of a cylinder embedded in R s and interact with a three dimensional SU(N) Yang Mills 
vector potential preserving a global chiral symmetry at finite N. As the circumference of 
the cylinder is varied from small to large, the chiral symmetry gets spontaneously broken 
in the infinite N limit at a typical bulk scale. Replacing three dimensional YM by four 
dimensional YM introduces non-trivial renormalization effects. 
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1. Introduction 

In the planar limit the number of Feynman diagrams grows only exponentially with or- 
der [|]]. At finite IR and UV cutoffs, one would therefore expect many observables to be 
given by a series in the 't Hooft coupling constant with a finite radius of convergence. We 
assume that this remains true in the continuum and thermodynamic limit. For asymp- 
totically free theories that would correspond to an expansion in scale [||. This has not 
been proven for pure SU(N) gauge theories, the topic of this paper. In both three and 
four Euclidean dimensions large distance phenomena are so different from short distance 
ones that one would guess that there are points of nonanalyticity in scale determining the 
radius of convergence of planar perturbation theory |J . If the theory has no qualitatively 
strong scale dependence, the radius of convergence might be governed by a singularity 
corresponding to some complex valued scale, an unphysical point. Our guess stems from 
the intuition that free-field short distance and long-distance confinement are qualitatively 



different to such an extent, that once we work in an approximation that produces a finite 
radius of convergence at short distances, a phase transition as a function of scale is plausi- 
ble. The simplest example would be a circular Wilson loop, expressed as a trace of a closed 
loop parallel transport in the fundamental representation, as a function of R the radius. 
For reasons we do not wish to be side-tracked into describing, it is unlikely that a critical 
radius exists in this case. However, if one looks at the distribution of the eigenvalues of the 
random unitary matrix whose trace is the Wilson loop a critical radius can be identified Q. 

Before one speaks about critical radii one needs to define renormalized eigenvalues, 
and that is somewhat complicated and involves the choice of order N/2 arbitrary renor- 
malization constants. A way to avoid the N dependence of the number of renormalization 
constants is to look at surfaces instead of loops and define an appropriate generalization of 
the Wilson loop operator || . One needs two (likely only one) new renormalization constant 
in 4D, and none in 3D. 

In this paper we study the 3D case and show that the new observable indeed has an 
unambiguously defined critical size associated with it. We then briefly look at the 4D case 
to identify in detail how it differs from the 3D case. Further work on 4D is relegated to 
the future. 

2. Setup 

We embed in R 3 an infinite cylinder S, defined by x^a),^ = 1,2,3 (a is short for a a ) 
with a = 1, 2. 



K LTTO\ K LTTG\ 

xi(a) = —cos— — ; x 2 (a) = — sin — — ; x 3 (cr) = a 2 ; —it < <j\ < it, -oo < a 2 < oo 

2tt l\ Itt L\ 



We define a two component gauge potential a a on the cylinder by 

dx 



a a = A li (x(<r))^± (2.2) 



The Dirac matrices on the surface are 



The new observable is : 

f\dA ]\d^dtb}e~^ f d3xTlF2+ ^ d2a ^^ [lad ° a ~ m ~ ilaaa ^^^ 
Q(m ' S) = " _i fd3xTlF 2 ( 2 - 4 ) 

The two dimensional massless Dirac operator is: 

D 2 (T,) = j a d aa - ilaa a {v) (2.5) 



1 2 + e, e > dimensional fermions coupled to four dimensional gauge fields were recently considered in 
a different context as an example of a theory with a conformal phase M . See M for related earlier work. 



Denoting the Hermitian generators of su(N) in the fundamental representation by T J , 
j = 1,..,N— 1 the nonabelian fermion current coupled to a is given by 



,P a (a) = ^(a) 7a T^(a) (2.6) 



The abelian vector current is 



J a = ^(a)^(a) (2.7) 

Power counting and symmetries allow for two local four- fermion counter-terms at m = 0: 

L\ = J^Ja, £2 = Ja.J a (2-8) 

Both terms can be simplified by a Fierz transformation. Both counter-terms are dimen- 
sionless, so there are no ultraviolet divergences worse than logarithmic. However, they 
won't be generated by the three dimensional gauge interaction. 

To identify possible ultraviolet divergences it is enough to replace the closed cylinder 
£ by an infinite two dimensional plane. The A^ propagator in Feynman gauge in Fourier 
space induces a two dimensional propagator: 
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(2.9) 



Taking / to vanish sufficiently fast at infinity, the limit A — > 00 leaves a finite expression, 

d2 P\\ ., 2, 1 



/ 
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7(PJi)-7= (2-10) 



unlike in 4D, where one is left with a logarithm of A. 

In short, because g 2 has dimensions of mass this mixed 2D/3D theory still is super- 
renormalizable. The main claim is that there exists a pure number, r c , such that at 
Af = 00 for Rg 2 N > r c we have 4(t^>) / while for Rg 2 N < r c , jj(4>ip) = 0. At 
finite N, -^{ipTp} = for all R. That -^(ipip) 7^ at infinite A^ means that the infinite N 
normalized single eigenvalue density of the random operator i.D2(£), p(/i), is nonzero at the 
origin: p(0) > 0. By clustering arguments the monomials ([jr'^ip] k ) can also be evaluated 
in field theory and expressed in terms of multi-eigenvalue densities of the operator Z?2(E). 

The upshot of these relations is that if p(0) > 0, the individual low lying eigenvalues of 

/ W\ 

iZ?2(S) are distributed in the same way as the eigenvalues of a matrix I , I with 

the probability density of the nxn matrix W given by exp(— ukTtW' W) where k is simply 
related to p(0) and n is taken to infinity ||. 

In this manner, the number r c identifies the smallest radius for which the smallest 
eigenvalues of — -D|(E) obey an explicitly known distribution law. The relation to random 
matrix theory does not depend on the chiral symmetry being non-anomalous |J. It only 
requires that the symmetry be obeyed by ^(E) gauge configuration by gauge configuration 
and that in the large N limit 00 > p(0) 7^ 0. It is also immaterial that at infinite N one 
can neglect the fermion determinant factor in the distribution of the gauge fields. But, 



this fact allows a more direct interpretation of Q(m, E) as a nonlocal operator in the pure 
gauge theory, similar to a Wilson loop, only associated with a surface rather than a curve. 
At infinite N the dependence of Q on m becomes nonanalytic in m at m = if g 2 NR > r c . 

3. An abelian exercise 

Replacing the gauge group SU(N) by 17(1) eliminates the gauge field nonlinearities and 
there is no parameter N to take to infinity. There are no phase transitions in the fermion 
world, but now the chiral symmetry is anomalous. As a result, for any R one expects a 
nonzero fermion condensate. In particular, we can replace the cylinder by an infinite plane. 
Chiral random matrix theory would be as applicable as before. In spite of the apparent 
non-linearity in the gauge field fermion coupling, the condensate can be computed by only 
carrying out Gaussian integrals. The calculation is almost identical to one way of solving 



the pure 2D Schwinger model [1C]. 

The fermions now live on the x% = plane and are not quenched. The condensate 
would be nonzero but also suffer from an infrared divergence if we quenched the fermions, 
just like in the Schwinger model |ll], |l2|| . 

The partition function is given by: 

Z= [dA„][dipd^}e ^o J "" J (3.1) 



Here a a = x a ,a = 1,2 are coordinates on the plane, e^ is the abelian gauge coupling. 
We wish to calculate 

(^(a)^(a')) (3.2) 

and obtain the condensate (ipijj) from the limit \a — a' \ —> oo assuming clustering. This 
strategy allowed us to set m = from the beginning. 

3.1 Decoupling the fermions 

The gauge potential on the plane can be decomposed as: 

A a = d a x + eapdpcf) (3.3) 

We now change Grassmann integration variables as follows: 

Here R, L denote the chiral components of the surface Euclidean Dirac fields. The change 
of variables exposes the anomaly and eliminates the gauge- fermion coupling. As a result, 
we can replace the fermion dependent part of the Lagrangian by 



S F = - d z a 



^flada^f + ^4>d 2 C 



(3.5) 



We now see that ipf,ipf are free and massless. The two dimensional scalar field <fi is 
determined by the electromagnetic field strength F\^ restricted to the u-plane: 

d a d a <t>(a)=F lj2 (a,x 3 = 0) (3.6) 



The fermion bilinears we are interested in are V, V denned by 

V = %b L ip R , V = $ R ip L , Tpip = V + V (3.7) 

In terms of the fields appearing in Sf we have 

V(a)V{0) = ${a)tf (0)^/ (<r)$f (0)e 2[ ^ ff) -^ (0)] (3.8) 

To compute the expectation value we need the portion of the eft dependent action coming 
from the pure gauge part. 

3.2 The gauge action for <fi 

If Sf = the expectation value of (f)(cr)(f)(0) is determined by the two point function of 
the field strength. The gauge action contribution to the (/> action is Gaussian with a kernel 
given by inverting the two dimensional kernel (Fi 2 (x)Fi 2 (0)) evaluated in pure gauge theory 
after setting X3 = and x a = a a . Using F^ u = d^A u — d u A^ we have the gauge invariant 
expression: 

(^12(^)^12(0)) = -^^3 + 5 -^- type of terms (3.9) 

47r \a\° a 

To perform the inversion we need the two dimensional Fourier transform of g(a) = 

1/|<t| 3 and require ultraviolet regularization. We choose to regulate by: 

(F 12 (a)F 12 (0)) = -| L l + ff{a) (3.10) 

47r (|cr| 2 + a 2 )2 2a 

Here a is a short distance cutoff, equivalent to 1/A where A was the high energy cutoff 
used earlier. The <5-function term is chosen to make the integral over a zero, because we 
expect to get a contribution to the expectation value of ( |3.8| ) only from gauge fields Fi 2 (a) 
which integrate to zero over the plane. 

g(k) = [ d 2 a w d ts = 2vr f°° ^^J {ka) 

J (|a| 2 + a 2 )§ Jo [a 2 + a 2 ]§ 

2nd [°° ado _ ., , 2ir _ ka 

J (ka) = — e ka (3.11) 



a da J y/a 2 + a 2 a 

Here, k = \k\. We now get: 

J d 2 a(F 12 (a)F 12 (0))e ll ° = || (l - e~ ka ) (3.12) 

We can take the limit of a — > at a fixed k and get 

J d 2 a(F 12 (a)F 12 (0))e^ = ^. (3.13) 

This result is compatible with (|2.10| ), obtained in Feynman gauge and using a sharp mo- 
mentum cutoff A at an intermediary stage, instead of the short distance cutoff a used here. 
Since 

F 12 (x) = -d 2 (j)(x), (3.14) 

the gauge contribution to the (ft action is 

Sgbfi) = \ I d2 ^( x ) i-d 2 } § <t>( x ) ( 3 - 15 ) 

e o J 



3.3 Condensate 



Combining the free fermion field propagators with the result of carrying out a Gaussian 
integral over <f> we obtain, 



(3.16) 



where 



Defining 



we obtain 
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2l 1 N 


\ f d 2 k l-e ik - a e§|cr| f°° 
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J J (2vr)2 kl + 2^ 47T J 



du 1 — Jn(u) 

2tt ^ " 



en C 



2vr ' 
2F(z) =7-ln2 + lnz + 



e -zt 

dt , = 

vTT^ 



(3.17) 

(3.18) 
(3.19) 



The dimensionless correlation function G{z) is given by: 



G(z) = ^(V(a)V(O)) 



1 



= 4F( 2 ) 



Jo 00 ^- 



47T 



, (27TZ) 

The integral in the above equation asymptotically goes like - and therefore 

2 

lim G(z) 



(3.20) 



47T 



(3.21) 



and we have a nonzero, finite, condensate. At short distances, the integral goes like —Inz 
and therefore the propagator goes like p-, the expected short distance behavior. 

No renormalization was needed and we got a finite condensate indicating that the 
spectrum of the two dimensional Dirac operator on the surface will exhibit chiral random 
matrix behavior as a result of its dependence on the background abelian gauge field. That 
gauge field is inherited from the 3D bulk, but its fluctuations are augmented by the feedback 
from the two dimensional fermions, via the anomaly term. 

4. Condensate in the non-abelian case 

We wish to show that one will get a nonzero finite condensate as above also in the SU(N) 
case and in the infinite ./V limit, so without the feedback of the fermion determinant. This 
can be done only numerically. 

We work on cubic lattices of L 3 sites with periodic boundary conditions on which live 
SU(N) gauge variables. The action is the standard single plaquette Wilson action and the 
lattice 't Hooft coupling is denoted by b. Large N reduction says that the infinite N limit 
for any L > L c (b) becomes L independent |p!q| . The L,b pairs at which we carried out 
simulations satisfy L > L c (b). 



We first put fermions on one of the 1 — 2 planes with antiperiodic boundary conditions. 



The overlap action couples them to the gauge link variables in the particular plane 16]. For 
each statistically independent gauge configurations, determined by the three dimensional 
Wilson action we calculate the two lowest positive eigenvalues Ai < A2 of the overlap oper- 
ator. We ascertain that there is chiral random matrix behavior and extract the condensate. 
This is repeated for several choices of b. We then check whether the condensate scales with 
b in the way expected of a physical quantity of dimension mass. We find that the answer 
is positive. 

The eigenvalues are computed using the Ritz algorithm and are used to estimate the 
condensate by formulas from chiral random matrix theory [13, 14 1: 
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(4.1) 

where U p is the parallel transporter around a plaquette and m w is the Wilson mass param- 
eter appearing in the overlap Dirac operator and V = L 2 . Chiral random matrix theory is 
diagnosed by the two determinations Xj agreeing with each other. 

Our results for the condensate are presented in Table ||. For b = 0.60 we have shown 
evidence for reduction since the results on 3 3 agree with 4 3 . Another example of reduction 
can be seen at b = 0.90 by comparing results on 5 3 with 6 3 . At b = 0.80, we obtain 
consistent estimates from A^ = 59 on 4 3 and N = 47 on 5 3 showing that we are in the large 
N and large volume limit. 
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0.59914(11) 
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0.68372(8) 


1.018(20) 
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47 
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1.020(15) 
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0.82256(1) 
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0.87594(0) 


1.024(12) 


0.0449(6) 
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1.50 


0.88458(0) 
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0.0437(5) 


0.0435(2) 



Table 1: Condensate for a plane, at several values of b and N on an L 3 lattice. 



Since the three dimensional coupling, b, 
has dimensions of length, we expect the con- 
densate to behave as ^ as one goes to the 
continuum limit of the three dimensional gauge 
theory if the 2D scale is set by the three 
dimensional bulk scale. The condensate is 
plotted as a function of the tadpole improved 
coupling in Fig. |l[ There is a curvature in 
the behavior indicating a sizable finite lat- 
tice spacing effect. To extrapolate to the 
continuum limit we used the so called tad- 
pole improved coupling 6/ instead of b (bj = 
eb). A fit to the sum of two terms proportional to bj and bj respectively agrees well 
with the data. 
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0.378(5) 


0.0681(10) 


0.0683(5) 



Table 2: Results showing the approach to chi- 
ral random matrix theory for a cylinder with a 
4x4 square base embedded on a 5 4 lattice at 
b = 0.7. 
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Figure 1: The data for the condensate is in Table [y. The b = 0.6 result on 3 3 lattice at N = 47, 
the b = 0.8 result on 4 3 and N = 59 and the b = 0.9 result on 5 3 and N = 47 have not been 
included in the plot. The fit includes a finite lattice spacing correction of the form b J . 



Two quarks on the plane will feel a linear potential at large distances reflecting con- 
finement in 3D SU(N) YM theory. We know what the condensate should be if this were 
a 2D SU(N) YM theory with the same string tension. When we compare the measured 
condensate to that of the hypothetical 2D theory we find that the condensate S in units of 
the string tension a, that is the ratio Y*/ ^fa, is 0.29 in the present case, whereas the exact 



value in purely 2D gauge theory is 0.23. In Casher's picture of spontaneous chiral symme- 
try breaking the condensate is determined by the size of a would-be massless bound state 



of two massless quarks in an s-wave [17]. A smaller bound state corresponds to a larger 
condensate. Within this picture we would then say that the size of these would-be bound 
state is smaller when the gauge forces have a 3D origin than when they have a 2D origin. 
In the 3D case the would-be bound state forms at a scale shorter than the scale at which 
confinement alone would have formed it. We can say this because in the 2D case the forces 
are purely confining. The comparison is possible because the kinematics of the fermions 
are identical in the two case. Since the 3D forces at shorter distances are strongly attrac- 
tive the fact that S/^/a is larger for 3D forces than for 2D forces makes sense. Actually, 
in this context, had we found the opposite relation, namely ( -5= ] < ( -§= ] , we could 
have claimed to have invalidated Casher's argument. Conversely, one could interpret the 
relative closeness of 0.29 to 0.23 as evidence that confinement makes a major contribution 
to the condensate in our 2D/3D model. 

5. The large N phase transition for the cylinder 

The lattice version of the cylinder has an s x s square as basis in the 1 — 2 plane. The 
fermions live on the surface of this tower and interact with the gauge fields on the links 
in the surface by the chirally invariant overlap action. The fermions obey antiperiodic 
boundary conditions round the s x s loops and in the 3 direction. 

For each statistically independent gauge configuration, determined by the three dimen- 
sional Wilson action we calculate the two lowest eigenvalues of the square of the overlap 
operator, < Af < A|, for the L 2 cylinders in one selected direction we call 3. The col- 
lected data produces histograms which are used to determine whether the distribution of 
the eigenvalues for given values of s, b becomes that of chiral random matrix theory in the 
large N limit, or, instead, indicates a positive spectral gap. In the former case, we extract 
a lattice value of the condensate as before. For a given s we shall find non-zero condensates 
so long as b < b c (s). To find b c (s) we fit the condensate to K^b c (s) — b in a regime of 
b- values which is close, but not too close to b c (s). We then check whether b c (s) scales as 
expected with s and find that it does, eventually producing a physical length for the critical 
length s of the base square of the cylinder. We looked at square bases of lengths s = 3, 4, 5. 

5.1 4x4 square base — Details 

Table || displays the approach to chiral random matrix theory at b = 0.7. The average of 
the eigenvalue ratio, (^), approaches the chiral random matrix theory prediction slowly 
and attains it only at N = 101. The estimate of the condensate obtained from the first 
eigenvalue changes little over the entire range of N values listed in Table || but the second 
eigenvalue takes much longer to converge. A plot of the distribution of the ratio of the 
eigenvalues in Fig. g exhibits a slow convergence to that predicted by chiral random matrix 
theory. 

We computed the condensate for several values of the coupling below b = 0.7 in order 
to study the approach to the critical point. Our estimates are listed in Table ||[ 




Figure 2: Distribution of j 1 at b 
embedded in a 5 3 lattice. 



0.7 and N = 13, 47, 101 for a cylinder with a 4 x 4 square base 
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0.0751(7) 
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101 
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1.00(1) 
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Table 3: Results showing the estimates for the condensate for a cylinder with a 4 x 4 square base 
embedded in a L 3, lattice and at various N. 



All results are consistent with chiral random matrix theory as the condensates predicted 
by the two eigenvalues agree with each other. A plot of the square of the condensate versus 
the bare coupling (we could have also used the tadpole improved coupling) in Fig. [3| shows 
that the behavior is linear indicating that the critical exponent is i. The critical value of 
the bare coupling for a 4 x 4 square base is estimated at 0.77(3). 

We studied the behavior of the two lowest eigenvalues at b = 0.80 for several values of 
N on a 5 3 lattice in order to see if it is in the symmetric phase. A plot of the data is shown 
in Fig. |j. The gap at infinite iV is estimated by a fit to a N~s (a standard soft edge random 
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Figure 3: A plot of the condensate listed in Table |] as a function of the bare coupling. 

hermitian matrix prediction) plus a subleading A r_1 correction which reproduces the data 
quite well as shown in Fig. Q. The fit forces the N = oo limits of the two eigenvalues to 
be identical. The estimated gap at infinite N is small, 0.00167(4), but still indicates that 
b = 0.80 is in the symmetric phase. The subleading term is essential for the fit to work 
and this indicates that b = 0.80 is close to the transition point. 

We now wish to compare the large N character of this transition to that we have studied 
for Wilson loops in the past. The essence of the latter was that close to the transition the 
eigenvalue separation in the least populated regime went as -k for large scales and as 



AT273 



for small scales, while, exactly at the critical size, it went as g /4 . Something similar 
happens to the level separation in our case: In figure || we show fits of the logarithm of the 
eigenvalue difference to a linear function of log(iV) with fitted slope and intercept. The 
fits cease being stable when subleading terms are added as the ranges of log(A2 — Ai) and 
log(iV) are not large enough. The two extreme couplings, b = 0.70 and b = 0.80 are close 
to the transition from different sides. Farther from the transition, in particular for smaller 
b, the slope in the fit would match better the appropriate random matrix expectation. As 
b goes through the transition the effective slope one gets from any fit to a set of data taken 
at finite N, would have to vary smoothly from -1 to -2/3. The change would be steepest 
at the critical b- value. One could define and effective b c (N) where N is the largest value 
in a well defined set used to get the linear fit and where the effective slope is exactly -3/4. 
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Figure 4: A combined fit of (Ai j2 ) at 6 = 0.8 and N = 13, 23, 29, 37, 47, 59, 79 for a cylinder with 
a 4 x 4 square base embedded in a 5 3 lattice. 



However, one should keep in mind that in the case of Wilson loops we dealt with 
an odd dimensional Dirac operator (living on the loop), while here we deal with an even 
dimensional Dirac operator. We expect that the correct random matrix models describing 
the low eigenvalues of the Dirac operator in the ungapped phase would differ in the two 
cases because chirality exists only in even dimensions. For odd dimensions, see JH 



5.2 Evidence for a continuum critical size 

In order to check whether the critical size also has a continuum limit, we ran simulations 
of cylinders with 3x3 and 5x5 bases. Our presentation of the 3x3 and 5x5 cases will 
be less detailed than that of the 4x4 case. 

We set N = 47 and used 4 3 lattices for the cylinders having a 3 x 3 base. Cylinders 
with a 5 x 5 base were embedded in a 5 3 lattice with N = 59, as the lattice spacing at 
the transition is smaller now. We computed the condensate for several couplings and the 
results are listed in Table HI All the data in this table are consistent with chiral random 
matrix theory. 

In order to see evidence for a continuum critical size we plot the square of the dimen- 
sionless condensate, £s, as a function of the dimensionless inverse size of the loop, bi/s, for 
an s x s loop in Fig. |6| for s = 3, 4, 5. The plot for s = 4 is the same as in Fig. || in terms of 



12 




InN 



Figure 5: Logarithm of lowest positive eigenvalues difference versus log(iV) and linear fits. The 
data is at couplings b — 0.70, 0.75, 0.80 for a cylinder with a 4 x 4 square base embedded in a 5 3 
lattice. 



the new variables. We see reasonable evidence for scaling and our combined estimate for 
the critical size is s c yfa = 1.37(10); for the string tension a, we have used the approximate 

formula y/abj = \ -^ IS]- 



6. Four dimensions 

For Wilson loops we know that the description of the large N transition is the same in 
2,3,4 dimensions. We have seen that the large N critical properties of the Wilson loop 
transition also extend to our new surface observable in 3 dimensions. So, we go to four 
dimensions, where renormalization is less trivial. 

6.1 Abelian exercise 

As before, we start with the an abelian exercise. The method of solution is identical to that 
used in the three dimensional case, so a brief description of the differences suffices. The 
difference comes in through the photon propagator restricted to the plane. Equation ( 3.12] ) 
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S X s 



1 



0.37674 \ A 2 



(IT 



Ei 



S, 



3x3 


0.45 


1.00(1; 


0.177(3) 


0.176(1) 


3x3 


0.46 


0.98(2; 


0.168(3) 


0.166(1) 


3x3 


0.48 


0.99(2; 


0.154(3) 


0.154(1) 


3x3 


0.50 


1.00(2; 


0.135(3) 


0.135(2) 


3x3 


0.52 


1.00(i; 


0.125(2) 


0.125(1) 


3x3 


0.54 


0.97(2; 


0.114(2) 


0.111(1) 


3x3 


0.55 


1.01(2; 


0.105(3) 


0.106(1) 


5x5 


0.67 


0.99(i; 


0.0954(14) 


0.0943(7) 


5x5 


0.70 


1.00(i; 


0.0881(13) 


0.0884(6) 


5x5 


0.72 


0.99(i; 


0.0822(12) 


0.0817(6) 


5x5 


0.75 


1.01(i; 


0.0785(12) 


0.0793(6) 


5x5 


0.78 


1.01(i; 


0.0680(10) 


0.0689(5) 


5x5 


0.80 


1.00(2; 


0.0650(10) 


0.0652(5) 



Table 4: Results showing the estimates for the condensate for cylinders with 3x3 and 5x5 
square bases embedded in 4 3 and 5 3 lattices with N = 47 and N — 59 respectively. 



is now replaced by 

J d 2 a(F 12 (a)F 12 (0))e^ 



e 2 Q k 2 \K l {ka) 



IT 



ka 



{ka) 



-^k 2 
Air 



bg(y) +0(ka) 



(6.1) 

The modified Bessel function K\(z), for z > is positive and monotonically decreasing. 
For z — > + we have 

K^z) = (z/2) log(z/2) + 1/z + O(z) (6.2) 



For z — > 00 we have 



K 1 {z) = s Jf z e-*[l + 0{z- 1 )] 



(6.3) 



The cutoff a cannot be eliminated. Keeping only the leading term in (6.1) we get the same 
expression we would have gotten had we used a sharp momentum cutoff in Feynman gauge, 
K = 2a- 1 : 



J V 



d 2 P\\d 2 p ± /(pjj 



pfi+Pi ' 



d 2 V\ 



< A 2 (2vr) 4 p 2 +p\ 7 p 2< A 2 
This impacts the correlation of the bilinears: 

1 1 



^L/(pjj)log[A/| P |||] (6.4) 



(V(a)V(O)) 



(2vr) 2 \a 



2 UI2' 



with 



F(\a\ 



1 f s du 1 - J (u) 



,4F(\a\) 



u 1 



; p = 2tt /e ; s = A\a\ 



(6.5) 



(6.6) 



log(u/s) 
We are interested in the large s behavior of F, now viewed as a function of s. 
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Figure 6: A plot of the square of the dimensionless condensate listed in Table y and Table as a 
function of the tadpole improved coupling. 



To this end we define: 

AF(s,p) = F(s)\ p=0 - F(s)\ p = P - f ^[1 - Jo(sx)] i 

Z Jq x p — 10 



log(x) 



Up to terms vanishing as s — > oo one has 



F( S )| p=0 = -[log( S /2)+ 7 ] 



and 



AF(s,p) = ^ log 



P 



log(s) 



6.2 Condensate on the plane 

We therefore find that 



A- 2 (V(a)V(0)) 



£. (A. 

4tt V2vr 2 



log(A|(7|) 



(6.7) 

(6.8) 
(6.9) 



+ terms that vanish as A|<r| — > oo 



(6.10) 

At finite A the condensate vanishes. However, the decay of the correlation function of bilin- 
ears with distance is slow. The coupling eo does not renormalize far from the fermion world, 
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where it cannot be screened. So, we cannot just eliminate the logarithmic dependence on 
A in ( |6.10 ) by introducing a A dependence in eo- 



To proceed in our search for a universal continuum limit we may investigate the effects 
of a Thirring coupling, the single four fermion term possible in the abelian case. 

6.3 Thirring term 

A Thirring term of the form 

g 2 

St = yjaja, ja = Tplatp, a = 1, 2 (6.11) 

is added to the action (the exponent of the integrand in the path integral is given by —St)- 
By a Hubbard-Stratonovich transformation this is implemented by shifting the elec- 
tromagnetic potential A a coupled to the fermion current to 

A a ^A a + B a (6.12) 

The action for A a is as before and that for B a is 

S B (B) = ^J dPaBl (6.13) 

B a lives only on the 2D plane. 
Writing 

B a = d aX ' + e a pdpct/ (6.14) 

we see that the scalar field entering the transformation to free fermion fields is now cp with 

V = 4> + ^' (6.15) 

Therefore the field entering the vertex operators V, V is p> and the transformation of 
integration variables from interacting to free fermion fields makes a contribution to the 
action given by 

S F = -1- [<pd 2 <p (6.16) 

For (ft we have the electromagnetic action we found before: 

* G = |Aog(iW (6 - 17) 

The Thirring interaction for (ft' is 

S T = -^ 1 J<ft'd 2 ( ft l (6.18) 

We now change the scalar field integration variables, (ft, (ft' to <p, (ft and integrate (ft out. 
The resulting interaction for ip is 



S M = 2 / > : 



d 2 d 2 1 d' 



n q 2 a 2 1 _ W , l 

eg iiog(-9VA 2 ) 



p (6.19) 
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Tracing this through amounts to replacing F in eq. (|6.5| ) by 

,-,/, in 1 f s du 1 — Jn(u) 

F(W\) = - — p ~ (6.20) 

■'O log(u/s)— g 2 p/ir 

where, s = A|a| as before. 

For p = we get the same expression as before, the Thirring term being shielded away. 
Our new AF is 

AF(s,p,g 2 ) = F(s)\ p=0 - F(s)\ p = P - f -[1 - J (sx)] \ (6.21) 

So long as g 2 is kept cutoff independent, up to terms vanishing at infinite UV cutoff, 
the answer is 



p 

AF(s,p) = -log 



■ log(s) 



_p(l+g 2 /ir) 
6.4 No continuum limit in the four dimensional Abelian case 



(6.22) 



We realize that both eo and g 2 do not renormalize and there is no physical infinite cutoff 
limit. The singular set-up of the problem in which the fermions are restricted to a zero 
thickness surface does not seem to have a universal description: the surface must be given 
some thickness and then continuum physics will depend on some of the details for how this 
was done. In three dimensions there was no such problem. 

In the four dimensional non-abelian case there is a new four fermion coupling which 
will renormalize and the gauge coupling is absorbed by dynamical scale generation in the 
bulk. We hope that these ingredients will provide for a nontrivial continuum limit in this 
case at the expense of one new free parameter. 

6.5 The necessity of a nonabelian Thirring interaction in four dimensions 

We now proceed to investigate the ultraviolet problem in the four dimensional nonabelian 
case, in the limit of infinite N. The fermions are restricted to the infinite 1 — 2 plane 
in four dimensional Euclidean space. The setup is a direct generalization of the three 
dimensional case. Working on an L 4 lattice we again find for all b < b c (L) that the lattice 
chiral symmetry is spontaneously broken and chiral random matrix theory agrees with the 
numerical data. Consequently, we can extract the condensate as a function of b. 

Our results for the chiral condensate are summarized in Table ||. 

We have quoted the logarithm of the condensate since we expect it to scale linearly 
with the coupling. Indeed a plot of the logarithm of the condensate as a function of the 
tadpole improved coupling, bj = eb, in Fig. |7| does show a linear behavior. The one loop 
beta function for the four dimensional theory would suggest a slope of — jj- and we obtain 
a value that is roughly half of this number. This is in sharp distinction with the situation 
in three dimensions. The exponential behavior might be explained by an effective Nambu 
Jona-Lasinio model. In four dimensions the condensate does not scale naively relative to 
bulk masses. The term C\ = JaJa from equation (|2.8| ) indeed is needed, as perturbative 
analysis indicates. With its introduction, a separate, adjustable, mechanism of dynamical 
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TV 


b 


e 


1 /Ai\ 

0.37674 \ A 2 ' 


In Si 


lnS 2 


13 


0.346 


0.52286(10) 


0.993(9) 


-1.893(10) 


-1.897(4) 


23 


0.346 


0.51513(5) 


1.009(7) 


-1.863(7) 


-1.855(3) 


29 


0.346 


0.51342(3) 


0.997(7) 


-1.842(7) 


-1.846(3) 


13 


0.348 


0.53097(10) 


0.996(9) 


-1.934(10) 


-1.938(4) 


23 


0.348 


0.52481(4) 


1.001(7) 


-1.905(7) 


-1.907(3) 


37 


0.348 


0.52290(2) 


1.000(7) 


-1.900(7) 


-1.900(3) 


13 


0.350 


0.53771(7) 


1.012(9) 


-1.985(10) 


-1.971(4) 


23 


0.350 


0.53239(3) 


1.004(9) 


-1.941(9) 


-1.941(4) 


13 


0.353 


0.54653(8) 


0.999(9) 


-2.016(10) 


-2.015(4) 


23 


0.353 


0.54183(2) 


1.013(6) 


-2.005(7) 


-1.992(3) 


13 


0.355 


0.55178(6) 


1.009(9) 


-2.045(10) 


-2.034(4) 


23 


0.355 


0.54741(3) 


1.005(9) 


-2.031(10) 


-2.026(4) 


13 


0.358 


0.55911(7) 


0.998(9) 


-2.074(10) 


-2.074(4) 


23 


0.358 


0.55508(2) 


0.999(7) 


-2.060(7) 


-2.062(3) 


13 


0.363 


0.57010(7) 


1.007(9) 


-2.160(10) 


-2.151(4) 


23 


0.363 


0.56646(2) 


1.004(7) 


-2.127(7) 


-2.127(3) 


29 


0.363 


0.56587(1) 


0.994(7) 


-2.121(7) 


-2.128(3) 


13 


0.367 


0.57810(6) 


1.009(9) 


-2.188(9) 


-2.172(4) 


23 


0.367 


0.57482(2) 


1.011(9) 


-2.143(9) 


-2.133(4) 


29 


0.367 


0.57412(2) 


1.015(6) 


-2.189(7) 


-2.173(3) 


37 


0.367 


0.57372(1) 


0.996(6) 


-2.179(7) 


-2.181(3) 



Table 5: Condensate for a plane computed at several values of b and N on 10 4 lattice. 

scale generation will become operative, restricted to the plane, and the scaling of the 
condensate will be affected by it. 

The numerical implementation of the C\ term with the required sign is a technical 
challenge. The renormalization of the 2D/4D mixed system will require substantially more 
work. 

We ignored taking the continuum limit for the time being and checked whether square 
cylinders in this case exhibited a large N phase transition of the same type as we have 
seen in three dimensions. Preliminary simulations indicate that this is the case. Only after 
gaining control over the continuum limit will we be able to address the question whether 
the structure of the large N phase transition is truly a feature of the continuum theory, as 
it is in three dimensions. 



7. Summary 

We have shown that Dirac fermions on a two dimensional surface embedded in a three 
dimensional SU(N) Yang Mills theory exhibit a large N phase transition separating short 
and large distance physics. The transition occurs on a cylinder whose base has a scale s 
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Figure 7: A plot of all the data for the condensate listed in Table g. A linear regression fit is 
shown as a solid line. 



when s is equal to a fixed number in units of the bulk correlation length. The exact value 
of the number will depend on the precise geometry of the cylinder base. 

We have also seen that the iV-dependence of the level spacing closest to zero goes from 



jr in the broken chiral symmetry phase to 



i 



AF273 



in the symmetric (gapped) phase. At 



the critical point, it goes as ~ J, 4 . In this sense the situation is similar to that found for 
Wilson loops. Also, the condensate vanishes at the critical size l c as ~ \/l — l c . Again, 
one can find an analogue of this in the Wilson loops case. 

We suspect similar facts are true in four dimensions, except that there one extra free 
parameter will have to be introduced in order to properly define the associated observable. 
Thus, the critical size l c of the base of a cylinder would depend on this extra parameter. 
We venture the guess that the exponents associated with / — l c and N will not depend on 
this extra parameter and will be the same as in 3D. Further work on the four dimensional 
case is needed. 
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